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Abstract 

The capture of trapped ultracold neutrons (UCNs) by closed hydrogenous impurities 
within a solid coating of the trap is discussed as a possible cause of observed anomalously 
large losses of UCNs in solid material UCN traps. Then significant losses of UCNs arise 
only if resonances occur in the UCN-impurity scattering amplitude. For a large size 
impurity, higher partial waves in the UCN-impurity interaction are important, and they 
are taken into account in the present paper. The method of the calculation is applicable 
to irregular shape impurities as well. A small distortion of an impurity shape, if it splits 
the resonance, can increase the UCN losses by a few times. UCN losses in the beryllium 
trap are calculated assuming they are due to the UCN capture by ice spherical impurities 
within the coating of the trap walls. Both s- and p-wave resonances contribute significantly 
to the UCN losses considered. As an example, observed anomalous large UCN losses are 
achieved if the average radius of the impurity is about 600 A and the impurity density 
is about 3 x 10^'^/cm^. A distortion of the spherical shape of the impurity could increase 
the UCN losses and therefore decrease the impurity density. 
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1 Introduction 



Losses of ultracold neutrons (UCNs) in material traps have attracted attention in recent years 
[1] . The matter arises for crystal materials such as beryllium and graphite [H El E] especially 
for beryllium [H [2] . For a liquid or solid Fomblin oil the observed losses could be explained by 
inelastic scattering processes [1] , but it can not [H [3], HI [5] explain the UCN losses in beryllium 
traps. Indeed, the UCN losses in beryllium traps depend slowly on the temperature T at low 
temperature [Il[2l[3]. By contrast, the inelastic scattering should lead to a ~ T'^ dependence 
of the UCN losses. The discussed losses cannot be a result of a surface heating caused by 
hydrogenous surface impurities pj because it requires too high a concentration of the hydrogen 
on the surface. In addition, in this case the losses should be significantly reduced by a degassing 
procedure, but the observed UCN ones are insensitive to it [Hill 13]. The discussed UCN losses 
cannot be from a coherent absorption of UCN by beryllium because the UCN capture cross 
section for beryllium is too small. The scattering of UCNs by vacuum cavities can increase the 
UCN losses ^ , but it can hardly explain the anomalously large losses observed [U [2] . 

In the present paper we consider UCN losses from the UCN capture by closed impurities 
within the solid coating of the trap walls. A low temperature of the trap is assumed so that 
only the elastic scattering and capture of UCNs can take place. Ice impurities are mainly 
kept in mind because the neutron-hydrogen capture cross section is large. Water within a 
closed cavity can hardly be removed by a degassing procedure. So losses are insensitive to the 
degassing process and do not depend on the trap temperature, nearly as it has been observed 
experimentally. 

A possible correlation of large losses of UCNs in beryllium with incoherent processes has 
been noted in lU [7] . In [1] the concept of a localization of UCNs around the lattice defect has 
been proposed which drasticly increases the UCN losses. This mechanism of UCN losses is not, 
however, acceptable since it is based on a mistakable solution of the scattering problem [8]. 

At low temperatures the losses of UCNs arise solely when absorption of UCNs takes place. 
If UCN losses are due to coherent absorption of UCN by the medium, then the well known 
expression [6] for the coefficient ^{y) of the UCN losses is as follows: 



where E is the UCN energy, U' is the real part of the optical potential, and rj = U" /U', where 
U" is the imaginary part of the optical potential of the medium. Inelastic scattering is negligible 
at low temperatures and so rj = b'/b, where b' is the UCN absorption amplitude and b is the 
UCN elastic scattering amplitude. Then for beryllium t] = 2.47 x 10"^ while the experimental 
losses [1] require 1] = 3 x 10~^. This exceeds the expected value of r/ by a factor of about 
100. The UCN losses due to UCN capture by a small size impurity are again represented by 
([T]) with the understanding that r] is replaced by pyw where pv is the portion of the total 
volume occupied by impurities, and w = W/U'. Here {—W) is the imaginary part of the UCN- 
impurity potential. For hydrogen w is known to be w ~ 3 x 10^^. Then rj ~ 10~^ requires 
j9y ~ 10 which is an obviously unacceptable result. 

The losses drastically grow for a large size impurity when a resonance in the UCN-impurity 
scattering amplitude occurs. Compared with the UCN losses by a vacuum cavity [B], they are 




(1) 



2 



increased by a factor, which roughly is waKU' j [Er]). In this case a is the radius of the impurity 
and K = ^J 2m{U — E)/h^ is the length of the UCN penetration into the coating; m is the 
neutron mass and % is the Plank constant (reduced). Under the discussed conditions this factor 
is ~ 100. 

As is shown in the paper, a relatively small impurity shape distortion, which splits the 
resonance into resonances, increases the UCN losses, roughly, by N times. In particular, a 
small distortion of the spherical impurity is able to increases the UCN losses from the /-wave 
resonance by about (2Z + 1) times in comparison with the UCN losses from a spherical impurity 
occupying the same volume. In the paper, however, as the first step, spherical shape impurities 
are mainly considered. 

Computations are performed under conditions where experimental data of the UCN losses 
are obtained [T]. Then the maximal UCN energy in the trap is about 23 cm, 38 cm, 46 cm, 52 
cm and 58 cm (this means that the UCN energy is the same as the gravitational energy of the 
neutron lifted to the given height; 1 cm energy corresponds to 1.025 x 10~^ eV). The radius 
of ice impurity is taken in the range (453.7 - 875.2) A. If the radius is less than 464.2 A, then 
a resonance does not occur for the UCN energy below 58 cm, and such size impurities do not 
contribute to the UCN losses. In the radius range being approximately (460-740) A, only an 
s-wave resonance occurs, and solely the s-wave UCN- impurity interaction is important. If the 
radius is larger than 740 A, then both s- and p-wave interactions contribute to the UCN losses. 
The s — p wave interference is taken into account, too. The UCN losses mainly arise from 
impurities lying rather far from the trap wall, but not from those lying on the trap surface. 
The d-wave interaction is negligible since a ci-wave resonance does not occur in the discussed 
energy-radius region. 

When an /-wave resonance is present, then the /-wave nonresonance background grows and 
contributes about (20-30)% to the UCN losses. Therefore a resonance approximation such as 
that employed in [6] is not used in the above computations. A large radius approximation [6] 
is not used in above computations, too. Nevertheless, both these approximations are discussed 
(in Sec. IV) for a semiquantitative consideration of the matter. 

An example of the distribution of the ice impurities is proposed fitting the experimental 
data [1]. Now experimental data are insufficient to fit them in the unique way. It assumes 
that no impurities have radii greater than 875 A. The impurity density is required to be about 
3 X 10^^/cm^. In fact the density falls when the size of the impurity increases. As a result, 
the average size of the impurity is about 600 A. The impurities occupy 0.25 of the total vol- 
ume. This is rather a large part of the volume, but not so much that the UCN capture by 
hydrogenous impurities should be discarded as a possible cause of the UCN losses. Perhaps, in 
accordance with aforesaid, distortions of the spherical impurity are able to decrease the volume 
occupied by impurities. At low temperatures the hydrogenous capture of UCNs remains to 
be a solely plausible cause for UCN losses, and so a knowledge of the cavity distribution in 
shape, orientation and size would be important. The method of this paper allows to perform 
calculations for more realistic impurity distributions. 

Under resonance conditions, the UCN-impurity scattering amplitude is able to achieve ex- 
tremely large magnitudes ~ 100km! In spite of this, the UCN can be treated as interacting 
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with a single isolated impurity. Indeed, a macroscopic effect from an interference of scattered 
waves each being formed by its own scatterer, could arise if many scattered waves interfered 
with the given scattered wave. Hence the effect could be mainly due to impurities separated 
from each other by a macroscopic scale distance ~ L in directions parallel to the boundary 
of the trap. Simultaneously, a sub-barrier scattered wave exponentially decreases going away 
from the scatterer. As the result, the discussed interference effect is negligible. Another effect 
could be that the scattered wave being formed by a given impurity, falls on another impurity 
after it was reflected from the trap boundary. The possibility of such a process decreases as 
independently of the magnitude of the scattering amplitude. So this process is negligible, 
too. This is demonstrated in Appendix A of the paper. 

The paper is organized as it follows. In Sec. II the wave function of the UCN in the trap is 
calculated. Inside the material of the trap it is represented as a superposition of wave functions 
ipQ. Each ipQ is the analytical continuation to sub-barrier UCN energies of the wave function 
of a neutron in the infinite homogeneous matter with impurities. By this construction, the 
wave function of the UCN in the trap is already matched on the boundary of the impurity. 
Matching on the trap boundary leads to a linear integral equation for the UCN wave function 
in the trap through ipo. The result is obtained without any restrictions on the UCN-impurity 
interaction. This is used in Appendix A to estimate interference effects discussed in the previous 
paragraph. They being negligible, ipo is given in terms of the relevant amplitude of the UCN 
scattering by a single isolated impurity in the infinite homogeneous matter, and so the UCN 
losses are calculated through this scattering amplitude, too. It is similar to the way in which 
the incoherent losses due to s-wave scattering of the UCN by the vacuum cavity have been 
calculated [6j in terms of the s-wave amplitude of the UCN-vacuum cavity scattering. 

In Sec. Ill the results of Sec. II are discussed for the spherical shape impurity case. 

In Sec. IV interactions of the UCN with a very small and a large size impurity are considered. 
In the large size impurity approximation, a semiquantitative consideration of the UCN losses 
is given, and a comparison with the case of the UCN scattering by a vacuum cavity [6], is 
performed. 

In Sec. V a small quadruple distortion shape impurity is discussed. 

In Sec. VI the UCN numerical computations are given, some details being given in Appendix 

B. 

2 Interaction of UCN with impurities 

In this Section a wave function of UCN in the material trap is calculated. Expressions for 
UCN loss cross sections and for coefficients of the UCN losses are given. It is implied that 
UCN interacts with material impurities, which presumably exist inside the coating of the trap 
walls. Basic results of this Section are applicable without restrictions on the UCN-impurity 
interaction. A set of impurity parameters is denoted as {a}. If UCN interacts with a spherical 
single impurity, then {a} = (rc,a), where is a radius- vector of the impurity center and a 
is the impurity radius. It is accepted that the plane ^ = is the border between the vacuum 
space {z < 0) and a material medium {z > 0) containing impurities. 
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The UCN wave function '?/'i(r, ko, {«}) in the vacuum space is given by 



V^i(r, ko, {a}) = e^^° ' + R{qo)e~'^'>^^"^'^'^°-^ + J dCq, ko, {a})e^^-^ e-''"'^'^^'d^q , (2) 

where r = { I ,z) = {x,y,z) is the radius-vector of the space point in question. In this 
case / is a two-dimensional (2D) vector in Xy-plane. Further, ko is the UCN wave vector: 
ko = ^o(q'o)) where go is a 2D vector in XF-plane. So go = fcosin6', where 6^ is a hide 
angle. Furthermore, -R(go) is the well known coefficient of the coherent reflection of UCN from 
the border, 

fco(g) +tK{q) 

In ([2]) and and throughout the paper, 

kl{q) = kl-q^, kl = 2mE/h^, Kl = 2mU/h^, [q] = + q^ ^ n'^ = ^1 ~ , (4) 

where m is the neutron mass, and U is the optic potential of the material medium; Rell > 
and Imll < 0. The sub-barrier case E < Rell is considered. The first term in the right side of 
(E]) describes the falling wave, the second term represents the wave formed due to the coherent 
reflection of UCN from the border, and the integral describes the incoherent wave due to the 
UCN-impurity interaction. The integration is performed from g = to g ^ oo. The condition 
q > ko determines the intrinsic reflection region. In this case fco(g) = i|A;o(g)|- 

The wave function ([2]) is matched at 2; = with the wave function '?/'2(r, ko, {a}) in the 
material medium off impurities as follows: 



^i(r, ko, {a}) = i)2ir, ko, {a}) 



dzi^iir, ko, {a}) = dzi)2{r, ko, {a}) 



2=0 



(5) 



2=0 



A full set of boundary conditions includes, in addition, the matching of wave functions on the 
impurity boundary. An important subtlety of the task is that V'2(r, ko, {a}) is represented by 
relatively a simple superposition of functions V'o(r, "p^, /t, {«}) each is an analytical continuation 
from {E — U) > Q io {E — U) <Q oi the wave function of the neutron in the material medium 
with impurities: 

V'olr, {a}) = e"^''" + ^o(r,"p, /€,{«}), b = ("p,m(p)), = -k^ ^ (6) 

where the flrst term in the right side describes the falling wave, and the second term arises due 
to the scattering of the neutron by impurities. In this case 'j) = {px,Py) is the 2D vector in 
XY-p\ane, p^ and Py being components of the UCN momentum, and k is deflned in Since 
the reflection from the border does not change the UCN energy, all ipoi^^ P > ^) {'^}) correspond 
to the same neutron energy E = U — hn^ /2m, but they are distinguished in "p^. Therefore 

V'2(r,ko, {a}) = y C2{'p,ko,{a})tpo{v,^,K,{a})(fp. (7) 
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If UCN interacts with a spherical, single impurity, then ipoi'^, ~P , i^i {c^}) ® is as follows: 

^ gi6|r-rc| 

^o(r, p", {a}) = e*^ '"^F(K, a, -ih ■ 



r — Fr 



where, as before, Fc is the radius-vector of the impurity center and a is the impurity radius. 
Furthermore, F{k, a, cos 6) is the scattering amplitude, 6 being scattering angle, cos 6 = h-h'/b'^. 
In this case b' is the final state wave vector [b = b' = in) which was replaced in ([6]) by its 
operator — iV. Outside the impurity the wave function ([6]) satisfies to the Schrodinger equation 
with the U interaction potential and has the true asymptotics at r — > oo as follows: 

tZ-oIf, ^, K, {a}) ^ e*''" + e**'"^F(K, a, bb'/fo^)^ 1 . (9) 

|r - Fc| 

The factor exp(zbFc) in the second term in the right side of ([6]) is the value of the plane wave in 
the center of the impurity. As it usually is, the scattering amplitude is expanded over partial 
waves amplitudes which are well known for a constant interaction potential [9J. 

Eq.([H]) is directly extended to the case of an arbitrary shape impurity. Then fq deter- 
mines a certain convenient point inside the impurity and F{k, a, cos 6) is replaced by a relevant 
scattering amplitude F(b , b', {a}). 

If V'oli") {c^}) is known, or if a reliable approximation to this can be proposed, then 
Ci{~q, ko, {a}) and C2{j) , kg, {«}) are calculated from boundary conditions For this pur- 
pose iI)q{y^~p ^K^{a]) in ([7]) is represented by means of Fourier integral over two-dimensional 
vector ~q as follows: 

V'o(r,"p,K, {a}) = e'^-^ i^o{'q,z,^,K,{a})d'^q, (10) 



where definitions are given in (jl]). Inasmuch as (fTOl) obeys the Schrodinger equation, then 
i^oCq, z, K, "p^, {a}) depends on z as exp {±'K,{q)z). Since ipoi^r, ~p , k, {a}) is a scattered wave 
formed in the right half space {z > 0), only exp {'K,{q)z) survives at z ^ 0. Hence 



dzM't, "p , K, {a}) = ^{q)Tp{'q, "p , k, {a}) 



:iii 



z=0 



where ilj{~q , J), {a}) = i^oCq, 0, J), k, {a}). Then ([5]) is turned out to be 

[1 + R{qo)]6\'q -q^) + CiCq , ko, {«}) = C2( p", ko, {a}) + B{~q , ko, {a}) , 

iA;o(go)[l - i?(go)]5^("g -%)- iko{q)Ci{~q ,^0, {a}) 

= -2{q)C2{'v, ko, {a}) + 2{q)B{-^ , ko, {a}) , (12) 



where 



5("g , ko, {a}) = j C2CP , ko, {a})i{~q , k, "p , {a}) (fp , (13) 
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and other definitions are given in (jl]), (fTOl) and (fTTl) . From (fT2l) . it follows that 



J?, ko, {a}) = 2fco(go) _ ^) _ , 

«;o(<?o) +««;(<?o) 

ko, {a}) = /'"^^^i ^ (Y, ko, {a}) , (14) 
fco(g) +m(g) 

where R{q) is given by The ~ (5^("g^ — go) term in C2("p^, ko, {a}) describes the wave 
penetrating into the matter from the vacuum, and the rest term is due to the reflection of the 
scattered wave from the boundary. Equation for B{ q , ko, {a}) is derived by substituting f|T^ 
to f|T3l) as follows: 

B{~q, ko, {a}) = [l + R{qo)]^i~q, %, k, {a}) - j R{p)B{~p , ko, {a})^("g , "p , k, {a}) d^p , 

(15) 

where definitions are given in ([3]) and (ITTl) . 

Once eq. lITS!) has been solved, the UCN wave function is calculated using ^\M . The total 
number A^;(ko) of neutrons lost per second is given by the integral of the UCN current density 
over the z = Q plane. Indeed, by first principals [9j, the change per second of the number of 
particles in a volume is given by the particle current going through the volume boundary. In 
terms of the wave function ([2]) the current density j( I ,ko, {«}) is given by 



j( / ,ko,{«}) = — Jm V^i(r, ko,{a})9^V^i(r, ko,{a}) 
m \ 



(16) 



where the right-top star denotes complex conjugation and the right part of (1161) is calculated 
at 2r = 0. Therefore, 

iV,(ko,{«}) = j j(T,ko,{«})rf2/ = hk,{q,)[{l - |i?(go)P)]5/m + iV,z(ko,{«}) , (17) 

where 5" is the area of the boundary, R{qo) is given by ([3]), and other definitions are given in 
(jl]). The ~ 5* term in the right side of f|T71) gives the number of UCNs lost per second due to 
coherent absorption of UCNs by the medium. The term A^j;(ko, {a}) gives the number of UCNs 
lost per second due to the UCN capture by impurities. Furthermore, Nuiko, {a})m/[/iA;o(go)] 
is the cross section crc(ko, {«}) of incoherent losses of UCNs in the trap. Indeed, by definition 
[9], the cross section of the process is the number of the questioned events per second being 
divided by the density of the falling flow. Using ([2]) and f|T6l) . one obtains c'"c(ko, {a}) as follows 

a,(ko, W) = 87r2/m[i?*(go)(-*)Ci(^,ko, W)]-47r2 / i^|Ci(^, ko, {a})pt/2g . (18) 

Jq<ko ko{qo) 

The last term in the right side of (ITSl) represents the cross section crs(ko{a}) of the UCN in- 
coherent scattering to the trap. Thus the first term is the total cross section of incoherent 
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processes 0. So (ITSl) is an optical theorem for incoherent processes under the discussed con- 
ditions. Eq. (fT8l) can be also written down using an amplitude As{~q ,ko, {a}) of the UCN 
incoherent scattering defined as follows: 

M^, ko, {a}) = -2mfco(g)Ci(^, ko, {a}) . (19) 

Taking into account that 

q = ko sin 9, dPq = kIcosO sin dOdcj) , (20) 
(7s(ko, {a}) can be represented as follows: 

a,(ko,{a})= / |A,(t,ko,{a})p.4^rffi., (21) 
J ko{qo) 

where integration is performed over the spatial angle Qg of the scattering of UCN into the trap 
{0 < 9 < vr/2). The factor ^0/^0(90) arises in ( l2Ti) because the velocity of the passing of UCN 
through the boundary is hko^q^j/m while the velocity of the returned neutron is hko/m. Then 
(IT8|1 is represented as follows: 

4:71 

(Te(ko, {«}) = ..—-Im[R*{qo)Mko, {a})] - a,(ko, {«}). (22) 
koiqo) 

The first term in the right side contains 47r/ko{qo) instead of the known An/ko since the velocity 
of the falling flow is hko{qo)/m. The R*{qo) factor takes into account a modification of the flow 
due to coherent interaction of the UCN with the matter medium. 

The losses Ni^ko, {a}) and Nu^ko, {a}) of UCN with given energy ko are obtained by the 
integration of the losses over directions of UCNs in the trap. Furthermore, Nii{ko, {a}) = 
4:7ikoac{ko, {a}), where ac{ko,{(y}) is the averaged cross section. It commonly is that an 
isotropic angular distribution of UCNs in the trap is assumed [U E] . In this case 



ac{ko,{a}) = — d(f)o a^ko, {a}) cosOosineo d9o = — crc(ko, {a}) — 2^ , 

4vrio Jo '^'^Jqo<ko k'o 



(23) 



where angles {6o,(j)o) specify direction of ko. 

Since interference effects from impurities are negligible (see Appendix A and Introduction), 
all the above relations are applicable to the interaction of UCN with a single impurity. For a 
spherical impurity the cross sections a^ko, {a}) and adko, W}) be denoted respectively 
as (Tc(ko, -2c, a) and ac(ko, -^c, «)• As before, Zc is a distance from the center of the impurity 
to the border 2; = 0, and a is the impurity radius. The macroscopic cross section is obtained 
multiplying crc(ko, Zc, a) by the impurity density which, generally, depends on a. Then it is 
useful in addition to introduce the macroscopic cross section for a certain convenient density 



^It is a certain mishmash in [6j where the cross section in the left side of is treated as the total cross 
section of incoherent processes, see eqs. (8) and (9) of appendix 6.13 in [6]. 
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no = This macroscopic cross section T,Qc(ko, Zc,cl) and the averaged cross section 

Soc(^o, Zc, a) are given by 

Eoc(ko, Zc, a) = no(Tc(ko, Zc, a) , uq = 10^^/ cm^ , Soc(fco, Zc, a) = rioadko, Zc, a) . (24) 

The coefficient of the UCN losses is the relation of the total number of the lost neutrons to the 
total number of neutrons falling on the boundary. For n = uq, the coefficient /io(ko,a) of the 
UCN losses, and the averaged coefficient Jlo{kQ,a) are given by 



yUo(ko,a) 



Soc(^05 Zc, a)dzc ■ 



(25) 



UCN losses arise only if the UCN interaction potential has the imaginary part somewhere. 
Indeed, the current density f|T6|) can be given in terms of the wave function \E'(r) in the medium. 
Off impurities \E'(r) coincides with ?/'2(r, ko, {a}) discussed above. This \E'(r) tends to zero at 
z ^ oo and obeys the Schrodinger equation as follows: 



V^*(r) = ^|K(r)-i?ol*(r) 



(26) 



where the ^(r) potential includes the impurity potential, too. As is usually done, multiplying 
both parts of fl26|) by \E'*( / , z), subtracting from the obtained equation its complex conjugate, 
and integrating the result over the z > half-space, one obtains that the total number NiCko) 
of neutrons lost per second is given by 



iV^ko) 







y ImV{v) 


^(r) 



(fr. 



(27) 



where j{ I ) is the current density at 2; = 0. If ImV{r) is zero all over (and therefore UCN 
absorption is lacking), then the UCN losses are absent. 



3 UCN losses from spherical impurity. 

The results of Sec. II are applied in this Section to the interaction of UCN with a spherical 
impurity. Hence ^o{r, "p*, k, {a}) is given by ([8]). Furthermore, 

e ""'I f ^■7i'.(i:-t)^-k(q)\z-Zc\ (^28) 



|r — rd J 271 K{q) 

Thus ip{~q , 'p ■, K, {a}) = "ipCQj 'p 1 0) in eq. lfTSl) is represented as follows: 
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where 



cos 'dip.q) = -li + K,{j>)'K{q)] I 



(30) 



and other definitions are given in (j4]). Furthermore, B{~q ^'k.Q, {a}) in (fT5l) is represented as 
follows: _^ 



B{ g ,ko,{«}) 



7iK{q){ko{qo) + m(go)) 



(31) 



where -FCg*, ko, Zc, a) is independent of Ic ■ Using (fT5|) . the equation for F(lj, ko, -Zc, a) is found 
to be 



F{~q, ko, 2;c, a) = F{k, a, cos ?9(go, q)) 



R{p) 



F(~p, ko, Zc, a)F{K, a, cos ?9(p, q))(fp . 



2ttk(p) 

By f fTSl) and flM|) . the macroscopic cross section of the UCN losses Eoc(ko,2c,a) is as follows 

''"(^°)^=F(^,ko,2c,a 



(32) 



Eoc(ko, 2;c, a) = 16?2ovr/co(go) 



Jm R*{qo 



q<ko 



fq 
27 



-2K{q)zc 



ko{q) + iK{q) 



ko{q) 



g-2K{go)^c^(^Y, ko, -2c, a) 



fco('?o) + i^^{qo) 



(33) 



As was discussed in the Introduction, in this paper UCN losses due to the UCN capture by 
impurities are mainly considered, ImU = being kept. Then (133|) is turned out as follows: 



Eoc(ko,Zc,a) = IGnon 



X 



/mF(g^,ko,Zc, a) - 2 



q<ko 27rKQ 



F{'q,ko,Zc,a) 



(34) 



Eq. fl32|) is easy solved when only the s-wave UCN-impurity interaction is taken into account. 
Then the UCN-impurity scattering amplitude is independent of the scattering angle, and so 
F( g , ko, Zc, a) is independent of the scattering angle, too. The solution of eq.(!32l) is easy found 
(cf. appendix 6.13 in |6j) being as follows: 



F{~q, ko, Zc, a) = Fq{k, Zc, a) 



1 



{K,a) - J{zc,k) 



J [Zc, k) = - 



-2k{p)zc 



Rip) 



k{p) 



pdp . 



(35) 



In addition, eq. (l32|) is analytically solved for the large size impurity when na » 1. Closed 
impurities being considered, then Zc/a > 1. Due to the exp[—2K,{p)zc] factor, only p ^ 
are significant in the integral in (l32l) . So F{'p', ko, Zc, a) can be replaced by F{0, ko, Zc, a). By 
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taking q = and calculating the integral over p at nzc >> 1, a simple equation for F(0, ko, Zc, a) 
is obtained. In this case F(0, ko, Zc, a) is as follows: 

F{0, ko, Zc, a) = e-^^-^cj^i^.^i 77 > (36) 

1 + ^-i^R[Q)F[k, a, -1) 

where -R(O) is given by ([3]) at g = 0, and a, — 1) is the backward scattering amplitude. 
Substituting (l36l) into (1321) . one easy derives F{~q, ko, Zc, a) for a general ~q, but it will not be 
used below. 

The foresaid equations of this Section are easy expended to the case of an arbitrary shape 
impurity. In this case a and Zc are replaced by a set {a} of parameters determining the shape, 
size and placement of the impurity and F{K,a,cosi!}{p,q)) is replaced by a relevant scattering 
amplitude. 

As is usually done for the spherically symmetrical potential, F{K,a,cosd) is represented 
through partial amplitudes Fi{K,a) as follows: 

l=oo 



r [K, a, cos 0) = ^(2/ + l)FKft:,a)PKcos^), (37) 
where P/(cos^^) is the Legender polynomial. In this case 

-iyprmp)/n)prmq)/n), (3^ 

where (pq is an azimuthal angle of ~q . For m > 0, x > 1, one finds that 



rn=—l 



Prix) = (v^^r-^ , PTi-x) = i-iyprix) , 

To solve fl32|) beyond approximations fl35|) and fl36l) . one represents ^(7?, ko, 2:c, a) as 

Z=oo m=i 



F(t,ko,^e,a) = 5^(2/ + !)(-!)' ^ e^^^^-'^-") ^ 7 go, Zc, a)Prmq)/^) • 

Z=0 m=-« ^ 

(40) 

From (|32|) . equation for Fi"^\zc, k, go, a) is found to be 



F["'\K,qo,z,,a) = Fi{K,a)Prmqo)/^^) - Fi{K,a) ^ C^r^^,, «:)F^^\«:, go, ^c, a) , (41) 



i'=|r?i 

where 



ClP'iz, K) = -{-lYC\-\z, -)y^T^^^{^1' + 1) , (42) 
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^11' ^) — ~ 



K{q) 



qR{q)Prmq)/K)Pirmq)/K)dq. 



(43) 



In this case R{q) is given by Q where fco(?) = ^|^o(q')| for Q > %■ Hence J{z, k) in fl5Sl) is none 
other than Cqq^(z, k). Both (P7|) and (155]) are below approximated by a finite number / < l^ax 
of terms @. 

For a constant potential, -F/(k, a) in fl371) is given by an analytical continuation in energy E 
to E = U — hn^ /2m of the corresponding partial amplitude [9] as follows: 



7r /«/iVi/2('««)'^m/2(H - ^^«+i/2('«a)-^/+i/2(H 

^l[K,a) = ( — 1) — — 

2 kKi+ii2{Ka)J[^^j^{ka) - KKl^^^^{Ka)Ji+i/2{ka) 



(44) 



where Ip{x), Kp{x) and Jp{x) are relevant Bessel functions; for any function /(x), it is defined 
f'{x) = df (x) / dx . Furthermore, 



= 2m{Ud + iW + E)/n^ = kI{(3 + iw) + kl = k'^ + iwnl; 
k'^ = Kl(3 + kl, i3 = Ud/U', w = W/U', U' = ReU, 



(45) 



where —{Ud + iW) is the potential of the impurity, and more definitions are given in In 
this case Ud > and W > 0. If ImU = 0, then real and imaginary parts of Clp\z, k) in (l43l) 
are given by 



ReC\-\z,K)^CZ\z,^) 



K{q)K^ 



, prmq)/n)ppmq)/K) 



-2 



^qe-'^^^''>J^^,Prmq)/K)Pn^q)/'^) 



ko '^0 



ImC[T\z, k) = Cnu/iz, K 



(m) 

211' l^? 



ko 



dq 



^0 



'e-'^^^''^^koiq)Prmq)/^)Pl^mq)/^) • (46) 



By using of ( l40l) . of ( 14TI) and of ( l46l) . eq.( l34l) can be transformed as follows: 

oo 

Y,[-i^Fr\^)]{-iy{2i + i) 



So.(ko,z.,a) = no^^^^e-^(->' 



1=0 



m=l 
m=—l 



r(/-m + l) 



r(/ + m + 1) 



Fi '{K,qo,Zc,a] 



(47) 



To derive (H7I) . one substitutes (HOj) to Then the second term in the right side of is 

represented as follows: 



oo oo oo 

' E E E 

m=—oo l=\m\ l'=\m\ 



6,, U'«^)(2^ + l)r(^ + ^+i) ^^''^^^^^ +l^r(/' + m + l) 



xFi^r\K,, go, 2;^, a)Fl;"''{K, go, ^c, a) • (48) 



(m) , 



^For this reason we do not discuss a convergence of the above series. 
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Using then eq. (HT!) and taking into account the first term in the right side of one obtains 
fH7|) . If fHT]) is cut off by / < Imax, then Fi^\k, Qq, Zf., a) is approximated as follows: 

^max 

fI'^\k, go, ^c, a)=Y^ z„ a)P,7(«;(go)/«) , (49) 

Z'=|m| 

where the ^;™'*( ) matrix obeys the equation: 

^["^(k, 2;, a) = Fi{n, a)6ii^ - Fi{k, a) ^ CP{z, k, a)A\7^^{z, k, a) , (50) 

l'=\m\ 

6jj' being the Kronecker symbol. Then 

poo Imax 

Mh,a) = no-rj / rf^; ^[-/mFf ^(«, a)](-l)'(2/ + 1) 

fTl—l { 1 I 1 \ ^'m.ax ^max 

X E tIt^ ^ 5: a)4?*(^, a) , (51) 

m=—l ^ ' ii = |rrt| /2 = |ni 

where c!2^j^{z, k) is defined in f H6|) . and the right-top star denotes complex conjugation. 

Throughout the paper, parameters in (jl]) and in (H5l) will be taken for the case of the UCN 
capture by ice impurities in the beryllium trap. As is commonly done, the optical potential 
of a neutron is 2Trh'^pnasc/'m where p„ is the nuclear density of matter and a^c is the UCN 
scattering length. The Be nuclear density is 1.235 x 10^^/cm^, and the ice molecular density is 
0.3175 X 10^^/cm^. The n — Be, n — O and n — p coherent scattering length is, respectively, 
7.79 fm, 2.9 fm and -1.87 fm. The imaginary part of the n — p scattering length is —4.63 x 10~^ 
fm. This is calculated from a neutron absorption cross section for the neutron velocity of 2200 
m/sec. The ice potential is the sum of the n — O and n — H2 potentials. The n — Be potential 
is found to be ?7 = 2.505 x 10~^ eV or in cm: Lq = 244.3 cm. The absorption of UCN by Be 
is neglected {Imll = 0). Then the parameters in (jll) and in f H5i) are found to be 

Ko = 1.1007 X lOVcm; P = Ud/U = 0.02452; w = W/U = 3.176 x 10"^ . (52) 



4 Small and large impurities 

To clarify main peculiarities of the UCN losses, the limiting cases na » 1 and na << 1 are 
discussed in this Section. 

If higher wave resonances (/ > 1) lie outside of the UCN energy range, then, as has been 
noted already and will be discussed below, eq.([25D is reasonable. When, in addition, absorption 
of UCN by the matter is negligible, then from flM|) the averaged coefficient Jlo{kQ,a) of UCN 
losses (1251) for the impurity density no is found to be 

~ fl \ A n J2{Zc,K,)ImFQ^{K) ~(o) 

Mh, a) = -47r no / ^ ^ r —da , J {Zc, k) = Qoo(^c, «) , (53) 
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where Fo{k, a) is given by ( l44l) at / = and C[^q{zc, is defined by (l46l) . In accordance with 
first principles ImFQ^{K, a) is negative. It follows from (l46l) that ImJ2{{zc, n)) is negative, 
too. Thus to calculate (l53l) in the leading approximation at IV ^ 0, only the leading term in 
ReFQ^{K,a) and in ImFQ^{K,a) are important. Thus ReFQ^{K,a) is given by fj44|) at Z = 
and ly = 0, while ImFQ^(K,a) is approximated as follows: 



IttiFq -^{k, a) 



w- 



K^(cot ka — ka/ sin^ A;a) 



2k [/tcosh(Ka) — cot fca sinh(Ka)]^ 
If the resonance does not occur, then rescattering can be neglected and therefore 



An 



d Zc J2(zc, K)/mFo(K, a) 



-^a 



If in addition K^a << 1, then from fl54p it follows that 

IniFQ^K^ a) = wo^Kq/S . 
Then, by using one obtains that 

2n 



ldoiko,a) 



riQwa 



Kq arcsin(A;o/Ko) — y i^o " 



(54) 



(55) 



(56) 



(57) 



Replacing no in (!57|) by the true impurity density and integrating the obtained expression over 
a, one comes the result multiplying by 4 to eq.([l]) with t] = pvw as was announced in the 
Introduction The UCN losses are small in this case. Much larger losses arise when there is 
a resonance in the UCN-impurity scattering amplitude as is shown below using the na » 1 
approximation. 

Only for simplicity one can assume in addition that fcoo >> 1 and k^a/ n >> 1. Then, by 
using ( l36i) in the ImU = approximation and after integration over go the averaged coefficient 
fl23l) of the UCN losses is found to be 



/Lio(/co,a) ^ 

Furthermore, d{e~'^'^'^'' /{2zc)) 
the integral over Zc gives: 



no 



dZr 



e ^'^'"'ImF^K, a, — 1) 



S-Zc 11 + 



22c 



-i?(0)F(K,a, 



(58) 



— (2^6 '^'^^'^ / {2zc)) dzc under the considered conditions. Thus 
47rno ImF{K,, a, —1) 



IJ,o{ko,a) 



21 n • cl'^*^'^'^) 

Ko/col-T (k, —a, 1) sm4| 

F(K,a,-l)i?(0) 



cos^ = Re- 



$(/t, a) = arctan 



/e~^^''\F{K,a,- 
V (2a)|sine| 



+ 



cos^ 
I sin^l 



\F{K,a,-l)\ 
cos^ 



— arctan ■ 



sin^l 



(59) 
(60) 



^Multiplying by 4 is made because the averaged cross section ([23| contains the 1/4 multiplication factor with 
respect to the average cross section defined in [U [6] . 
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The backward amplitude F{k, a, —1) is calculated by (137|) through partial amplitudes (jHj), each 
being approximated at >> 1, as follows: 



.F,(.,a)^(-l)'e''''-''^-l^''''b-^'l^-"'l^\ C,^JI±1I^, (61) 

K + kcot{ka-7rl/2-ai/2) 

where /^//ta < 1, ka » 1 and na » 1. Keeping ai is reasonable only if / >> 1. The 
resonance is determined by the condition 

K + kcot{ka-7Tl/2-ai/2) = 0, (62) 

where k is k at W = 0, cf. ( H5i) . Therefore, nearby the resonance where a\kl — k'^ \ / k << 1 , the 
backward amplitude F{k, a, —1) is given by 

4(2/ + l)A:,^(l+/3)e^--"e-"' 

where kj = 2mEr/h'^, = Kq — fc^, and Er is the resonance energy. The background is denoted 
by ellipsis. Other definitions are given in (jl]) and fH5|) . Then the resonance contribution to the 
coefficient of the losses (|25|) is found to be 

/^o(/co, a) ^ , , . „ ^ ,,,, . , = , (64) 



kr\sm^\^{e,-k^^y + 



WK, 



where $(k, a) is calculated by fl60p with the understanding that the whole F{K,a, —1) is sub- 
stituted into (1^ . its background being included too. If na — > oo, then mainly the P ~ Ka 
terms contribute to the background. So F{n,a,—1) ~ aexp(2Ka). Furthermore, the imagi- 
nary part of the partial amplitude fl6T|) is given by its derivative with respect to W. Hence 
ImF{K,a, —l)/\ReF{K,a,—l)\ ~ {wnla/k), where w and k are defined in fH5l) . So the reso- 
nance losses (El]) prevail the background ones when o^I/cq — << 1. It follows from ( l62i) that 
under the discussed conditions many resonances arise, each being approximately determined by 
the condition ka ~ 27rs provided 2Tis/{ka) « 1, and where s is an integer including zero. At 
given / the distance between neighboring resonances is 2TTk/a which is much greater than l/a^. 
Resonances for / and / + 4n where n is an integer number, are close to each other as ~ l/a^. 
Hence the resonances do not overlap in the region a^\kQ — k'l\ « 1, and so the total resonance 
contribution to Jloi^o^ cl) is a sum over resonances. 

The discussed asymptotics is not fully achieved for a < 880 A considered below. In this 
case a single s-wave resonance occurs in the UCN-impurity scattering amplitude, and a single 
p-wave resonance is added if a > 740 A. Then a nonresonance piece of exp(— 2Ka)F(fi;, a, — 1) 
is rather ~ than ~ a. Hence in the a^l^Q — >> 1 range both terms on the right side 
of (1601) cancel each other and therefore $(/t, a) decreases. In this case rescatterings cease to be 
important. The resonance losses flMl) prevail the background ones when a\kQ — k'^\ << k^. 

In spite of the fact that rescatterings from the wall, generally, spread the resonance, the 
width of the resonance (!64|) remains ~ wHq. This is due to the fact that the main contribution 
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UCN reduced energy UCN reduced energy 



Figure 1: Coefficient Ho{ko,a) of UCN losses at qkq = 6.5 against the UCN reduced energy 

to fl58|l arises from so great Zc that T2{zc) ~ max[|fcQ — fc^|,wKg]. Indeed, approximating the 
integrand in fl58|) by its resonance piece, one obtains that 



~ , r 47rrir2(;^c; , 



C ; 



2^ 



2 8(2/ + l)fc3(l + /j)/,e^--("---) 
r, = ^Ko r,(.J = 2z.a.g(.g + /5A:,^) " ^''^ 

The expression under the integral is non other than the resonance piece of the UCN capture 
cross section adko, Zc, a) introduced in Sec. II. For simplicity, / ~ 1 is considered. The integral 
diverges exponentially up to so great Zc that T2{zc) ~ max[|/i;Q — /c^|,wkq]. Hence these Zc give 
the main contribution to the integral as was announced above. 

The shape of the resonance peak assigned to (|64l) is quite different from the shape of the 
Breit-Wigner resonance. Indeed, contains 1 / ^/{^'^^^k^y~+JvJK^ , but not 1/[(/cq — + 
[wKq)"^]. Besides, $(fi;,a) is significantly altered when k^ runs within the resonance interval. In 
particular, a) is asymmetric with respect to the top of the resonance. The main features 
of Jlo{ko, a) are demonstrated by Fig.l where Jlo{ko, a) for ukq = 6.5 is shown against the UCN 
reduced energy y = E/U = k^/nl. Its detailed behaviour in the resonance range is shown, too 
(the right-side figure). 

UCN losses per second are determined by the integral over /cq of the expression which is 
Jlo{kQ,a) multiplied by the magnitude of the UCN flux J^(fco), see Sec. VI for more details. 
The leading contribution to the integral gives the range: wkq << |/cq — kj\ « Since 
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(66) 



This expression can be also derived directly from fl65l) . Indeed, the integration of fl65l) over fcg 
gives the result: 167r^rir2(2;c)/(r2(zc) + Ti)- This expression is integrated over Zc by means of 
introducing r2(^c) as the integration variable, and eq.( l66i) arises. 

The relation of the resonance losses to the background ones roughly is ln[l/?i;(/toa)^] which 
is ~ 10 under the conditions considered. Generally, the resonance in the partial amplitude 
(16T|) increases an average nonresonance magnitude of ImFi{K,,a) over the integration interval 
that increases the background piece of the UCN losses. If the resonance does not occur, then 
the partial amplitude (!6T1) has a zero (in the W = limit) in the integration interval. It, 
generally, reduces ImFi{n,a) that decreases the losses. As a result, partial waves of the last 
type contribute to the losses about a percent or smaller. These rough estimations are confirmed 
by the numerical results of Sec. VI. An accuracy of fl59|) has been estimated to be about (20- 



It is instructive to compare fl66|) with the losses from the scattering of UCN by a vacuum 
cavity [B]. By contrast to foresaid, in the last case Ud = W = in fH51) , but the imaginary 
part U" = 7]U' of the U potential is taken into account. Examining ([3SD with Ud = W = and 
at OK — i> cx), one concludes that the resonance piece of the UCN losses is again given by (!65|) 
with the understanding that Fi is Fi = r]k^/{aKr) and, in addition, that F2(2;c) and e{zc) are 
calculated by (16T1) at /5 = 0. For Z = it agrees with the expression (17) of appendix 6.13 in 
[6] taking into account the difference in the definition of the average cross section and up to 



5 Possible increasing of UCN losses for non-spherical 
impurities 

The leading contribution to the UCN losses at large a is given by the sum over resonances 
of fractional losses fl66l) . If the resonance is degenerated in some quantum number, then it 
is represented by a single term in the sum. Indeed, the numerator and the denominator in 
(|6^ each is proportional to the number of the degenerated resonances. In particular, each a 
resonance in the /-wave partial amplitude is {21 + l)-degenerated with respect to the azimuthal 
quantum number, but its contribution flM|) to Jlo{ko, a) contains no {21 + l)-multiplier. A 

^ In [6], there is a mistakable extra factor 2 in Fi and certain insignificant inaccuracies in e and 
^The losses explicitly depend on the cavity radius a as 1/a, but not a'^ seemingly occurring in eq.(20) of 
appendix 6.13 in [6] 



30)% 




are now as follows: 



(67) 
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spherical shape distortion is able to split the resonance into the {21 + 1) resonances. Then the 
UCN losses due to this resonance could increase by about (2/ + 1) times. Nevertheless, the 
volume occupied by impurity could remain about the same. It is demonstrated below with 
an example of a small quadruple distortion of the spherical impurity when the shape of the 
impurity is given by 

r = a + ^EpsnpUs, ess = 0, Sps = e^p^ns = Vs/r, = 1, el^ « a? , (68) 

where the Eps set gives the shape distortion. The summation over twice-repeated indices is 
performed. In this case is the s-component of the radius-vector r, and the center of the 
spherical impurity is placed in the point of origin. For the sake of simplicity, only mixing of 
states within the / = 1 multiplet is considered. The wave vector of the falling wave is kept to 
be (0,0, m). If an^ » 1, then the / = 1 piece 'ipmatij) of the UCN wave function outside the 
impurity, and its derivative V'mat(r) with respect to r, are given by 

^™at(r) ^ -^ns + fn— -, V'LtW ~ -^^^3 - fn^— . (69) 



Only the exponentially large term is kept in the falling wave in fl69|) for the calculation of 
the resonance piece of the scattering amplitude fn in the leading approximation. The I = 1 
piece ipimpi''^) of the UCN wave function inside the impurity nearby the boundary fl68l) . and its 
derivative ip'^j^p^r) with respect to r are given by 

, , cosA;r , ,/ / x fcsinfcr , , 
V'imp (r) ^ ^—dn, ^imp{r)^ dn, (70) 

where d = {di ,d2 , d^) is an independent of r. 

At Ers = 0, the resonance magnitude kr of the wave vector is found from eq.( l62l) at / = 1 
(and «! = 0). To find the leading term in the backward amplitude, the linear in Ers terms on 
the boundary fpSl) need to be kept only in (1701) . Then equations matching the wave functions 
(169|) and (!70l) on the impurity boundary (168!) . follows (r = 1, 2, 3): 

Se'"^ e"''" cos ka k sin ka 

— — 0,^3 = Jr H dr Ersds , 

Zna a a oa 

Se*^" e~'^"'K k sin ka k"^ cos ka 

-T, "r3 = Jr H dr H z £rs"s ; (71) 

Za a a oa 

where 6rs is the Kronecker symbol, and the summation over twice-repeated indices is implied. 
From (jTTI) . the backward amplitude F{k, a, —1) ~ — /s about the resonance (!62l) is found to be 

^, 12fc,^(l + /j)e^--"det33[i^,. + 2fc,%,/(5a)] ~ ^ , 2 _ , 2 .7,^ 
^' ' aiK^r+(3k^r)det[ix + twKl)6,s + 2k^Ejs/M]' ° ' ^ ' 

where det[xSjs + 2krEjs/5] is the determinant of the matrix whose elements are given in the 
square brackets (j, s = 1,2,3) of the above expression, and det^slxS js + 2krEjs/5] is (33)- 
minor of the determinant. To obtain the coefficient of the UCN losses, the amplitude fl72p is 
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substituted into (159!) . One can see that (172|) has the resonance, if x = Xj = —2k'^ej/{5a) where 
Ej is an eigenvalue of the {ejs} matrix. If a^\xj — Xs\ >> 1 for any discussed Xj and Xg, then 
the UCN losses will be the sum over the losses from every resonance, they will be three time 
more than the UCN losses in the e^s = case. The volume Vim occupied by impurity increases 
as follows: 

Vim = J d COS ed(f) Y (^o + ^ersrirnsj ^ -vra^ V ^ ^ 8^) ' ^^^"^ 

Hence the relative change of the impurity volume goes to zero when a — >■ oo. So, in the a ^ oo 
limit, the distortion of the impurity shape can increase the UCN losses remaining the impurity 
volume being about the same. As of now, it has not been studied whether the discussed 
increasing of the UCN losses takes place in a realistic range of a. 



6 Probability of UCN losses 

In this Section, there is given UCN losses from ice, spherical impurities calculated under con- 
ditions where the losses have been measured experimentally. An example of a radius impurity 
distribution is proposed which fits the experimental losses of UCNs in the beryllium trap. As it 
was noted in the Introduction, now experimental data are insufficient to fit them in the unique 
way. The proposed fitting must be considered only as an preliminary example. The calcula- 



ho 


58 


52 


46 


38 


23 


l{ho) 


0.4867 


0.4608 


0.4334 


0.3940 


0.3065 


ar{l, ho) 


5.109 


5.424 


5.791 


6.391 


8.120 


ar{2, ho) 


7.758 


8.206 


8.728 


9.583 


12.05 



Table 1: Reduced wave vector 7(/io) and ar{l, ho) impurity radius. 

tions are performed for a narrow cylindrical beryllium trap [Ij, its radius R being R = 38 cm, 
and its length L being L = 14 cm. Numerical data (1521) are imployed. 

For each of five discharges pQ the height ho of the trap is respectively 58 cm, 52 cm, 46 cm, 
38 cm and 23 cm. Simultaneously, ho measures i?max in centimeters. The reduced wave vector 
kmax/i^o assigned to -Emax, will be denoted as 7(/io) where ho is measured by centimeters. As 
an example, 7(58) is the reduced wave vector assigned to ho = 58 cm. The reduced impurity 
radius ano is kept within the range 4.994 < qko < 9.633 where the low limit corresponds to 
the s-resonance being at 7 = 7(58) + 0.01 = 0.4967, and the top limit corresponds to the 
d-resonance at 7 = 7(63) + 0.01 = 0.5173. An /-wave resonance is occurs at E = -Emax at a 
certain reduced radius of the impurity. This reduced radius will be denoted as ar{l, ho) where 
ho is measured by centimeters. The reduced radii in the region of interest are given in Tab.l. 

Experimentalists [Ij determine a probability of the losses per second r~^(i?max) of UCNs 
with energies up to the given maximal energy -Emax- In line with the foregone text, it is useful 
to introduce in addition TQ^{E^i,x,a), which is probability of UCN losses per second caused 
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Impurity reduced radius Impurity reduced radius 



Figure 2: Probability TQ"^(_E'max, o) of UCN losses against the impurity reduced radius in the 
region where resonances are not available; the trap high is 52 cm (curve 1), 46 cm (curve 2), 
38 cm (curve 3) and 23 cm (curve 4). 

by the UCN capture by impurities of the given radius a and of the density uq = 10^^/cm^. 
Hence ro"^(-Emax, a) and r~^(i?max) are given by 

1 f hk 

To\Erai,^, a) = —— / p{ko,h) — -Jlo{E,a) kodkodS , (74) 

^"^(^max) = / no— — ro~^(£'max, o)c?a, (75) 
Jo da 

where Jlo{E,aKo) is an averaged coefficient of the UCN losses fl2Sl) . nodn{a) is the impurity 
density in the (a, a + da) range, and p{{ko, h) is the energy-space density of UCNs in the trap 
which will be discussed below. In this case h is the height h of UCN over the base of the trap 
The integrations in ([71]) are performed over the surface of the trap and over k^ from fco = 



to fco = kmax = \J '^rnE^^^/h^ . Furthermore, At^ Ni{E^e.x) is the total number of UCNs in the 
trap as follows: 

A^i(^max) = / p(A;o, h) kodko dV , (76) 



where the integrations are performed over the volume of the trap, and over ko from /cq = to 

A common assumption [H [6] is that inside the trap UCNs have the isotropic velocity dis- 
tribution. Furthermore, due to the Earth gravitation field, the UCN energy is related with the 
height h of UCN over the base of the trap as mgh = -Emax ~ E, where g is the acceleration of 
the gravity. So the UCN density p{ko, h) in both the space and the momentum space is given 
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impurity reduced radius impurity reduced radius 

Figure 3: UCN losses rQ"^(_Emax, a) versus the impurity reduced radius ano for 4.99 < an^ < 7.6 
where only s-resonance is available and for 7.6 < ano < 9.6 where s- and p-resonances are 
available; the trap height is 58 cm (curves 1), 52 cm (curves 2), 46 cm (curves 3), 38 cm (curves 
4) and 23 cm (curves 5). 
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s-wave only. 
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4.953 


4.834 
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(s+p)-waves 


4.826 


4.844 


5.031 


5.136 


5.165 


5.140 


5.077 


4.992 


4.907 


4.940 



Table 2: Comparison of UCN losses/sec. Tq ^(i^max, x 10^ at ii^max = 58cm due to s and 
(s + p) (the last line) interaction in the region where p-resonance is unavailable. 

by a (5-function type expression as follows [H [6] : 

p(kQ,h) = ciS{mgh - E^s,y, + E) ko, (77) 

where ci does not depend on the UCN wave vector and on the UCN location. An explicit 
expression of rQ"^(-Emax, o.) for the considered trap is given in Appendix B, see eq. flB.2l) . 

Fig.2 shows UCN losses in a region where resonances do not occur. As it has been noted 
already, in this case the losses are very small, ro~^(-Emax; o,) ~ (10~^ — 10~'^)/sec. If only a s-wave 
resonance occurs in the scattering amplitude, then solely the s-wave UCN-impurity interaction 
is important. It is demonstrated by Tab. 2 at -Emax = 58 cm where the losses due to the s- 
wave UCN-impurity interaction are compared with the losses calculated when, in addition, the 
p-wave UCN-impurity interaction is taken into account. 

If a resonance is present in the UCN-impurity scattering amplitude, then UCN losses increase 
as it is demonstrated by Fig.3. When the resonance firstly occurs (it is an s-wave resonance) at 
E = E'max, a high peak arises. Then the losses fall, but remain rather large, r~^(i?max) ~ (0.5 
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Figure 4: UCN losses Tq ^(-Emax, o.) versus ukq in the region of the picks; the trap high is 58 cm 
(curves 1), 52 cm (curves 2), 46 cm (curves 3), 38 cm (curves 4) and 23 cm (curves 5). 



- 1) X 10~'^/sec. The losses sharply increase again when the second resonance (it is an p-wave 
resonance) occurs. Then the losses fall, but remain to be ~ (1 - 2) x 10~'^/sec. In Fig. 3 UCN 
loss probabilities are also separately presented for off peak regions (the bottom figures). Peak 
tops are not seen in Fig. 3 because of the peaks are extremely narrow. The UCN losses within 
peaks are shown in Fig. 4. The UCN losses are quite large in the peak ranges, but the integral 
contribution from the peak to the UCN losses is not prevailing because of the extremely narrow 
width of the peak. 

Experimentally measured UCN losses p[j can be fitted, for instance, by means of the impurity 
piecewise-smooth radius distribution as shown in Fig.5. In this case an impurity nondimensional 
density dn{a) / {dano) is given against the impurity reduced radius a/to- The impurity density 
in the {a,a + da) range is Uq Kq da dn{a) / (dano) . As above, tiq = 10"^^/cm^. Eqs. fl74|) . fl75l) and 
(IB.ip are employed in the calculation. In Tab. 3 UCN loss probabilities T~^(£'jnax) calculated 
for the impurity radius distribution in Fig.5, are compared with the experimental data [IJ. The 
impurity density n within the radius range considered is calculated as follows: 



n = uq 



dn{a) 
da 



da = 2.85 X lO'Vcm^ , 



(7^ 
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Figure 5: An example of the impurity nondimensional density dn{a)/{daKo) fitting Serebrov's 
data. 
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23 


Serebrov's data 


21.8 


17.4 


15 


11 
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Theor 


20.7 


17.9 


14.2 


11.4 


7.9 



Table 3: Serebrov's data (the second line) compared with UCN losses under an impurity radius 
distribution given in Fig.5 (Theor). 

where amin — 4.994/ko and a^ax — 9.633/ko- The average radius Uav corresponding to Fig.5, 
is found to be aav = 6.33/ko = 575 A, and a portion py of the volume occupied by impurities 
is pv = 0.25. As it was discussed in Sec. V, a distortion of the spherical shape impurity is 
potentially able to reduces pv- 
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Appendix A Interference effects from impurities 

Effects from the interference of scattered waves are briefly discussed here. The s-wave USN- 
impurity interaction is only taken into account. 
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Impurities being in number, ipoi''^, 1^, {(^}) & is given by 



M^,^,{a}) = y2Ai- (A.l) 

i=l 

where ri is the radius-vector of the i-th impurity center, and the Ai set is calculated by the 
matching of the wave function off impurities with the UCN wave function inside each an im- 
purity as follows: 

e'P-^--''^P>-F^ + F^Y. sA,e-^'-/rr,, = A„ , n = 1 . . . , iV ; r„„ = |r„ - , (A.2) 

where Fn is the UCN scattering amplitude on the n-th single impurity. So the A = {Ag} 
column is found to be 

A = {F-'-y)-'u, (A.3) 
where u = {«„}, F = {F^m} and y = {ynm} are given by 

"-n ; ynn ^ nm '-'nm^ m ; 

ynm = e~'^''""'/rnm for u =^ ui; ymm = . (A.4) 

Thus the effect of n-th and of m-th impurity on each other is negligible when Fnynm and 

^ ^ ^ 

FmUnm. both are small. Even when F„ = 100 km, Fnynm is less than 0.01 already for |/„ — lm\ = 

1.5 X 10^^ cm. The macroscopic effect is due to the interference from those scatterers, which 

separated from each other by a macroscopic scale distances being much more than the distance 

above. So An can be approximated by the relevant amplitude of the UCN scattering on the 

single isolated impurity. 

To discuss the interference under the reflection of the UCN from the trap boundary, one 

~ > 

notes that in the discussed case, B{q, ko, {a}) is found from eq. (fT5l) with the understanding 
that now the function ipoi'^y "p"? {o^}) is given by flA.ip . In this case B{q, /cq, {«}) is found to be 

B{q, ko, {a}) = V V ^ Bs.' , (A.5) 

s=i s'=i 7rK{q){ko{qo) + m(go)) 

where Bgg' = Bg's and the matrix B = {Bgs'} satisfies an equation as follows: 

B = {F-^ - x)-^[l - (Z'^) + J)B] . (A.6) 

Hence B is found to be 

B = {F-^-x-J^''^ -J)-\ (A.7) 
where matrix elements of matrices J = {Jmn} and J*^'^-' = {Jmn} are given by 
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where C^{zn,K) is given by ( l43l) . Due to a singularity at = in the integrand, Jmn 
decreases at | — /^l — > oo, as follows: 

^ pf^O I I 

2(1 + t)koe-^^^'-+'"^ ^ ^ , (A.9) 

-/ |2 

that is nonexponentially. To obtain (\A.9h one first integrates in J^n over the azimuth angle 0, 
keeping | Z„ — /ml — * cxo. The obtained integral is represented as 



where ellipses denote the integrand in the last integral. The above integrand is the same as in 
the first term. The last integral decreases exponentially when | — Z^l ~^ oo. The calculation 
of the first integral leads to ( lA.lOl) . 

By using ( JTSl) . ( |23ll and (1A.5I) . the averaged cross section adko, {a}) of the UCN losses in 
the y = approximation is found to be 

?c(fco, {«}) = E (^^-^^'^ + ^^^)- 1^"' - '^^'^ - Jl^'j^^n' n V ' ^^■'^^^ 

m,n=l s=l 

where an averaging over each / ^ is performed, as well. Thus the A correction in adko, {a}) 
due to the two impurity interference is as follows: 



^ m,n=l 



47r ImA^ ImF^ VmnP + ImJ^n ImF^ ^ 2/?e[Jm„(F„ ^ - Jn ) TT ^^^^ ^ a i o^ 



s=l 



From flA.gp . a relative correction flA.121) to the leading term of adko, {o;}) is roughly found to 
be ~ N/{SkI) ~ ndo/^o, where do is the length of the trap coating, do ~ 5000 A, and n is the 
impurity density. This correction is extremely small for any reasonable n. 

So, (lA.lip is a sum of the losses over the UCN losses from each a single, isolated impurity, 
as it is considered throughout this paper. 



Appendix B UCN losses in a cylindrical trap 

To obtain ro"^(-E'max, in the case of interest, cylindrical coordinates (p, 0, z) are employed. In 
this case 2;-axis goes along the length of the cylinder laying horizontally, p is a minimal distance 
from the given space point to 2;-axis, and ip is an angle in the perpendicular to z-axis plane, 
-0 = at the bottommost point of the trap. The infinitesimal element dS of the side surface is 
dS = Rdzdip, and for each of the butt-end, dS = p dp dip. The integration over ip is performed 
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employing 5-function in (!77j) . The integration with respect to z over the side and with respect 
to p over the butt-ends are performed without difficulties. Then (1771) is found to be 



2U 



nN{ho)LKo 



J ^0 



2RL 



Ll^RVLl - [kl/Kl - ho/Lo + R/Lo] 



where 



(B.l) 

LQ = U/mg , ho = E^s.x/mg . (B.2) 

In the calculation of A^i(emax) one integrates over E employing the (5-function in (177|) and 
integrates over z. In doing so dV = pdpipdz. The result is as it follows: 



N{ho) 



Lq V Lq 



— T "■V'xl [~ cosip — — costpo) , 



(B.3) 



where 

R-RcostlJo = ho. (B.4) 

The integration in fIB.SP is performed keeping the radicant being positive and in addition, 
ip < ipo and p < R. The integral can be calculated through Legendre function Pylil — h^/R) 
as follows: 



N{h, 



R\ 



5/2 



- 



vr 



LoJ V2 



1-1- 



ho 
R 



2-\ 



p,';,{i-ho/R) 



(B.5) 



Below eq.( ]B.5[) is proved for ipo < 0. The calculation for ipo > is performed in the same 
manner. 

For ipQ < 0, there are two integration regions in and p being as follows: 



<\i)\ < liJol , < p < R (regionl) and 
\ipo\ < IV'I < ^ ; < p| cosT/^l < R\ cost/'oI (reg.2). 



:B.6) 



The integral over each a region will be denoted respectively as Ni{ho) / L'J'^ and N2{ho) / L^^"^ . 
The integration over p is performed using that the indefinite integral 



p dp\/ {p cosip — -Rcost/'o) 



is equal to 
Hence 



2p 



dip 



3 cosip 
3 cosip 



(p cos ip — R cos ipi 



N3/2 



[cosip — cosipoY^'^ 



15 cos^ Ip 



{p cosip — R COS ipii 



N5/2 



(B.7) 



:b.8) 



15 cos^ Ip 



COS Ip — COS IpQ 



l5/2 _ 



cosV^o)^^^ 



(B.9) 
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and 

N,{ho) = 4 / #- — (- cosV^o)'/' • (B.IO) 

J 15 cos^ tp 

The second term in flB.91) is integrated by part using that dtp/ cos^ ip = dtsmip. Hence 

Niiho) = / dip cosip{cosip - cosipof^'^ sm'ijjo{- cosipof^"^ . (B.U) 

3 J 15 

Furthermore, 

2i?5/2 r i?5/2 



(i-?/' cos?/' (cos 'i/' — cos'?/'o)^''^ = TT ^ sin^ ■?/'o-Pi/2('^Q^'^o) ' (B-12) 
where is the Legendre function. The N2{ho) integral is easy calculated, and as a final 



result, f B.5 ) arises. 
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